Energy absorption of quasineutral plasmas through electronic edge-modes 
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Ultracold quasineutral plasmas generated in the laboratory are generically inhomogeneous and ex- 
hibit small charge imbalances. As will be demonstrated, via a hydrodynamic theory as well as 
microscopic simulations, the latter lead to efScient energy absorption at the plasma boundary. This 
proposed "edge-mode" is shown to provide a unified explanation for observed absorption spectra 
measured in different experiments. Understanding the response of the electronic plasma compo- 
nent to weak external driving is essential since it grants experimental access to the density and 
temperature of ultracold plasmas. 

PACS numbers: 52.35Dm,52.55.Dy,52.27.Cm 



In recent years ultracold plasmas, produced by photoion- 
ization of laser-cooled atoms [T], Q , have opened up a new 
realm of avenues to study various plasma phenomena, 
ranging from collective behavior 0, Q and strong corre- 
lations to low-energy atomic processes These 
studies in the ultracold regime rely on experimental ac- 
cess to the plasma parameters, such as the density and 
temperature of the different plasma constituents. 

With optical techniques such as absorption imaging Q 
or fluorescence measurements one can directly image 
the positions and velocities of the ions in the course of 
the plasma expansion. Tracking the dynamics of the elec- 
trons has turned out to be more subtle due to their lack 
of internal structure. One central approach is to use a 
weak radio frequency (rf) field and record the yield of 
electrons which escape due to resonant rf-heating. One 
can deduce the density [HI and temperature [lH of 
the electron plasma, provided the relation between the 
plasma parameters and the resonant frequency is known. 
In [lo| this was originally done by assuming resonant ab- 
sorption at positions r with density n(r) for which the 
field frequency lo matches the local plasma frequency 



a;p(r) = ATTe^n{T) 



(1) 



Based on this model the average absorp tion wa s con- 



V 47re^n, solely 
(n2)/(n). The 



eluded to be peak at a frequency Wp 
determined by the average density n 
density evolution extracted in this way, was found to be in 
excellent agreement with theoretical calculations [l^ 
for a wide range of initial plasma parameters. Subse- 
quently, it was suggested [l5| that the plasma has a con- 
tinuous excitation spectrum and absorbs energy through 
a damped quasi-mode of frequency ~ O.STcDp. While un- 
derestimating the densities found in [loj this approach 
gave a good description of subsequent rf-absorption mea- 
One may conclude that neither the pic- 
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surements 

ture of local absorption nor the quasi-mode permits a 
consistent interpretation of existing experiments, which 
would facilitate a reliable probe of ultracold plasmas. 

In this letter, we resolve this seeming discrepancy and 
reveal the existence of an additional mode that provides 
a unifying consistent description of all previous experi- 
ments. Our results are found to be consistent with the 



general characterization of eigcnfrequencies in an inho- 
mogeneous plasma [l6| . To make contact with experi- 
ments and to account for coUisional damping, we also 
present a hydrodynamical description that confirms the 
existence of the derived mode and illustrates the role of 
the quasi-mode for energy absorption. Finally, to avoid 
any phenomcnological parameters, we have performed 
large-scale molecular dynamics (MD) simulations with 
up to 50000 electrons. All three approaches give results 
in excellent agreement with each other, providing an ex- 
haustive characterization of the identified eigenmode and 
its role for energy absorption. 

The existence of the eigenmode is based on one prop- 
erty which has not been taken into account so far, namely 
the fact that the experimentally realized plasmas are 
generically only quasi-neutral. They exhibit a charge im- 
balance 



(5 = (TV, - N,)IN, 



(2) 



between the number iV^ of ions and of electrons. This 
imbalance results from a prompt loss of electrons upon 
creation of the plasma and, more importantly, from elec- 
trons that continue to leave at later times due to the pres- 
ence of small electric fields during plasma expansion [17 1 . 
The charge imbalance creates an outer edge of the elec- 
tron density, which is, consequently, described by a trun- 
cated Gaussian distribution. As we will demonstrate, 
this leads to dominant energy absorption at the plasma 
edge and produces an absorption-resonance frequency in- 
timately linked to the charge imbalance, lo = l^{5). 

The "edge-mode" , as we will call it, can explain the 
energy absorption reported in experiments [id . [T]| and 
may also be suitable to reconfirm the existence of Tonks- 
Dattner like modes at finite temperature [ll|, whose 
properties in unconfined plasmas are yet to be established 
firmly. 

We start with the theoretical description of a charged 
particle density at zero temperature, n(r) = ne(r) + 
(5n(r)e~*'^*, where ne(r) = no/(r) is the equilibrium den- 
sity and 5n is the perturbation induced by driving the 
plasma with a weak external rf field, Ee~"^*. Substitut- 
ing this density into the underlying hydrodynamic equa- 
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tions yields to leading oder in Sn 

,2 



UJ^ + UjJV 



/(r) V20-V/(r)V0 = -EV/(r) , (3) 



where Wpo = '^p(O) is the peak plasma frequency and the 
rate v models the effects of coUisional damping. The den- 
sity perturbation has been replaced by its corresponding 
electrostatic potential (/) through Sn = V^0. Without 
driving and damping (E = 0, = 0), Eq. ([3|) reduces to 



V(eV0) = 0. 



(4) 



which is simply Maxwell's equation for a medium with 
dielectric function 



e(r) 



1 



W2 



/(r) 



1 



^p(r)^ 



(5) 



In the following we focus on the spherically symmetric 
case and assume that the equilibrium density rte(r) is 
bound to within a sphere of radius R. This implies that 
the radial potential 4)i{r) represents a perturbation of 
well defined angular momentum I which has to fulfill the 
boundary condition 



e(r) 



dine 



dr 



dine 



dr 



(6) 



r=R 



at the interface R, with (pf denoting the potentials for 
r>R. 

For the case of a homogeneous density fir) ^ Q{R — 
r) , Eq. (U) reduces to a "free" Laplacian with general 
solution (f>e{r) = cir^ + C2r~^~^ . Since must be finite 
at r = and for r — > oo, we obtain <j)f = cir^ and 
= C2r^^^^ . Inserting these solutions into Eq. © gives 
w/wp = {(■/{(■+ l)]^^'^- These are the familiar multipole 
excitation modes for a homogeneous spherical plasma, 
where the dipole mode (£ ~ 1) has uj = a;p/\/3. 

For a general form /(r) Eq. Q is solved numerically 
by propagating the logarithmic derivative of (j^f outwards 
from very small r where the "free" solution (f)^ oc r^ 
dominates. At the boundary R of the plasma the prop- 
agated solution must be matched to the free outer so- 



lution 



which yields the desired eigenmode 



frequency w. 

In Fig. [T] the dipole excitation mode {£ = 1) according 
to Eq. (jl]) is shown as a function of the charge imbalance 
S, which we can relate to the cutoff radius R. Since 
the plasma electrons tend to neutralize the Gaussian 
shape ion cloud from the center, the resulting electron 
density distribution is described by a truncated Gaus- 

sian distribution = e^2^Q{R — r), where the ra- 

dius R is determined by the number of electrons through 
Ne ~ 47rn(0) / drr^fji{r). As one can see, the modes 
form a smooth function between the two limits of a neu- 
tral and a highly imbalanced plasma. The latter case 
((5 — > 1) implies a nearly homogenous plasma within a 



small sphere of radius i? -> for which the eigenmode is 
the dipolar plasma frequency as derived above and indi- 
cated by a horizontal line in Fig.[TJ 

In order to clarify the role of the derived edge mode for 
energy absorption we include a finite damping rate v and 
consider the corresponding energy absorption density per 
unit time 



psir,Lu)^^Rc{F3*) 



\F\' u;pir)^u 



(7) 



Here J is the current density and F = E — Vc/) the total 
electric field. To obtain ps{r,uj) for a given 6 one must 
solve Eq. ([3]) which can be done similarly as described 
before: The logarithmic derivative of is propagated 
outwards well beyond the plasma edge R where it can 



be matched to the free solution 



cx r 



To propagate 



across the plasma edge we have softened the so far ide- 
alized sharp edge by introducing a distribution function 
for r > i? which connects smoothly with /fl(r) at R, 



/;j(i?)(e('-«)M-l) + l ' 



(8) 



where the length A controls the softening of the edge. 

In Fig. [2] we show the resulting spatially resolved ab- 
sorption spectrum for a neutral plasma (S = 0) and for a 
plasma with a finite charge imbalance. As can be seen, 
local energy absorption prevails in both cases. Integrat- 
ing ps{r, uj) over r gives the energy absorption rate Pa (w), 
which is shown in the left part of the figures in Fig. [2] For 
6 = (Fig. [5^) this resembles the absorption spectrum of 
[Tsl , which peaks at the quasi-mode frequency. Note that 
this picture provides a unifying frame for the two models 
put forward in [lo| and [I3 , with the only difference be- 
ing the strength of local energy absorption at w = Wp(r). 
Along this line, one immediately understands the finite- 
d modification of the quasi-mode absorption, which can 
be estimated from the neutral-plasma absorption spec- 
trum po{r, uj) via the truncated integration (dashed lines 
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FIG. 1: (color online) Excitation mode of an ultracold plasma 
for different charge imbalances S (Eq. ^) from solving the 
eigenvalue problem Eq. Q. The dashed line is the limit of a 
homogeneous spherical plasma, the dashed-dotted line indi- 
cates the local plasma frequency Eq. ([T]). 
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FIG. 2: (color online) Double differential power density 
log[4:-Kr^ps{r,uj)] for a neutral (5 = 0, (a)) and a charged 
plasma {S = 0.2, (b)). The damping rate is = 0.007 and the 
edge has been smoothed over a range A/cr = lO""* in (b), see 
text. 



in Fig.ED 



R{5) 



(9) 



However, at finite S there appears an additional sharp 
maximum extending above the quasi-mode background 
(cf. Fig.[2b). This additional resonance originates from 
strong energy absorption at the plasma edge where lo- 
cal neutrality is maximally violated. Therefore, we call 
this resonance in a spatially truncated plasma an "edge- 
mode". In fact, the edge- mode is exactly the eigenmode 
CO derived above from Maxwell's equation for = and 
v = as can be seen in Fig.|31 

To further substantiate these findings and to deter- 
mine appropriate values for the phenomenological damp- 
ing rate v, we have also performed molecular dynam- 
ics simulations of a large electron plasma confined by 
the Gaussian neutralizing ion background. As seen in 
Fig. [3] the simulations confirm all our findings outlined 
above. Moreover, the enormous peaks (note the loga- 
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FIG. 3: (color online) Energy absorption of an ultracold 
plasma externally driven by a weak rf field of frequency lj 
for different charge imbalances 5. Vertical lines represent the 
eigenmodes from Eq. Q, the solid lines Pi(t<j)/max[Po('^)] 
are a numerical solution of Eq. (|3]) with v = 0.007, and the 
circles connected by lines represent AE{iu,5)/ ma,x[AE{(jj,0)] 
from an MD simulation with 25000 atoms. 
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FIG. 4: (color online) (a) Ion background (dashed) and elec- 
tron density for a neutral {S = 0) and a charged plasma 
{S = 0.2) from an MD simulations of 50000 atoms at three 
different Debye screening length, Xd/ct; (b) The influence of 
screening on the edge-mode. 



rithmic scale in Fig. [3]) suggest that an ultracold plasma 
with a significant charge imbalance absorbs dominantly 
through the edge-mode, even at non-zero temperature. 

We can estimate the effect of finite temperature by 
identifying the smoothing A (cf. Eq. (|5])) with the Debye 
screening length Ad ~ y^ksT/ {e'^Ann). This is demon- 
strated in Fig. 2^ with our molecular dynamics simula- 
tion, where we show several density profiles for different 
temperatures, expressed through Xjj /a. Indeed, one can 
see how the plasma edge broadens according to the inter- 
play between 6 and Ad/ct. As a consequence, the edge- 
mode absorption peak gets broadened and weakened 
with increasing temperature, i.e., \d (Fig.13}^). However, 

— 1/3 

for typical experimental values Xd/o' ^ 2N^ <^ 1 
and therefore edge-mode absorption still dominates 
quasi-mode absorption. 

In the experiment [lol [Tl| energy absorption is mea- 
sured by recording an enhanced evaporation of resonantly 
heated electrons, which accelerated towards a detector 
by a small extraction field i?oxtr during the course of the 
plasma expansion. As established in independent opti- 
cal measurements [Tsj , the density indeed stays Gaussian 
during the plasma expansion, which proceeds in a self- 
similar fashion according to (T{t) ~ Vco + v'^t'^, where 
the expansion velocity = kBT{0)/AI is given by the 
initial electron temperature T(0) and the mass M of the 
ions [l, [ll] . Since the density consequently decreases in 
time, the maximum electron yield appears at different 
times with changing rf-frequency and indicates resonant 
absorption at different densities. However, since the ex- 
traction field i?cxtr constantly extracts electrons during 
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FIG. 5: (color online) Edge mode (red) of an ultracold plasma 
as a function of charge imbalance in comparison with experi- 
mental results [l^l (□) and [ll[ (x), see text. The quasimode 
response {S = 0, ^'Hd its remnant for finite S is also 

indicated (black). 

plasma expansion, the charge imbalance also changes dy- 
namically. Hence, in order to determine max[P(a;, 6)] as 
predicted by our analysis (see Figs. [T] and [3]), we need to 
know S{t). This task is tremendously facilitated using the 
result of a recent investigation on the influence of the ex- 
traction field on the electron escape dynamics [TtI , where 
it was found that the electron yield Ncsc{t) is a univer- 
sal function of one parameter a = E cy^t,- a {tY / (An Ni) for 
times t > tovap after initial evaporation. This allows us to 
determine S{t) = Ncsc(t)/Ni by optimizing a single quan- 
tity, the extraction field -Eextr, which was only estimated 
in the experiments [l^, [HI ■ The result of this analysis 
is shown in Fig. [5] Most importantly, both experiments 
consistently fall on the same universal curve which co- 
incides with the edge-mode for increasing S. For small 
charge imbalances (5^-0 the experimental data seems 
to approach the quasi-mode limit. Note, however, that 
data points at very small S are missing since they are 
collected at times close to the initial evaporation where 



it is not possible to determine S{t) without additional ex- 
perimental information. This uncertainty, however, does 
not change the present conclusion or impair the use of 
our theory to unambiguously probe the electron plasma 
through dynamical measurements of the absorption res- 
onance and the charge imbalance during the expansion 
of an ultracold plasma. 

To summarize, we have presented an exhaustive the- 
oretical analysis of ultracold plasmas with a trun- 
cated non-homogeneous electron distribution and demon- 
strated perfect agreement with large-scale molecular dy- 
namics simulations. The truncation resembles charge im- 
balances, a generic but so far scarcely discussed feature of 
realistic ultracold plasmas. Using a macroscopic hydro- 
dynamic formulation we have identified the energy ab- 
sorption at zero temperature as a resonant edge-mode at 
well defined eigenfrequency in contrast to the continuous 
quasimode characteristic of a non-homogenous plasma 
without edge. Our analysis demonstrates that the edge- 
mode absorption takes place at the plasma boundary, 
which supports resonant electron escape as observed ex- 
perimentally. The identified absorption mechanisms pro- 
vide a unified description of existing experiments, open- 
ing the way for a precise probe of ultracold electron plas- 
mas by rf-hcating measurements. 

The notion of the truncated electron distribution could 
also explain additional resonances found experimentally 
above the fundamental mode as true "Tonks-Dattner" 
resonances in the ultracold electron plasma since the lat- 
ter is actually confined by the charge imbalance. How- 
ever, substantiating this conjecture requires in future 
work a careful analysis of temperature effects on the 
plasma edge which after all provides a dynamical and 
not a static confinement. 

We thank Dan Dubin, Kevin Twedt and Steve Rolston 
for enlightening discussions. 
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